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Abstract. - We investigate the influence of nonstationary 1//^ noise, produced by interacting 
defects, on a quantum two-level system. Adopting a simple phenomenological model for this 
noise we describe exactly the corresponding dephasing in various regimes. The nonstationarity 
and pronounced non-Gaussian features of this noise induce new anomalous dephasing scenarii. 
Beyond a history-dependent critical coupling strength the dephasing time exhibits a strong 
dependence on the age of the noise and the decay of coherence is not exponential. 



The question of the phase coherence of a quantum two level system (TLS) in a noisy 
environment has motivated numerous recent works. In particular it is of crucial interest in 
the context of quantum computing with solid state quantum bits [1]. In this case it is now 
believed that dephasing of the TLS by a low-frequency (1//) noise is a limiting mechanism 
which certainly deserves further studies. The microscopic origin of these low-frequency fluctu- 
ations depends on the considered system and is mostly not understood [2] . Hence, in previous 
works on this problem the 1// environment was modeled phenomenologically either by a set 
of harmonic modes {e.g. spin-Boson model) [3] or by an ensemble of independent bistable 
fluctuators (Dutta-Horn model) [4,5]. The first model assumes a Gaussian distribution for 
the fluctuations and stationarity of the noise. As this low frequency noise appears to be 
an essentially out of equilibrium phenomenon, this last hypothesis is questionable [6]. The 
Dutta-Horn model does not suffer from this restriction [5] , but assumes the presence of a broad 
distribution of relaxation times for the independent fluctuators, or equivalently a flat distri- 
bution for their random energy barriers. However, both the nature of these fluctuators and 
the origin of these long relaxation times are beyond the current experimental knowledge. We 
consider the interactions between the fluctuators as a possible source for these slow dynamics. 
, In this letter, we explore the consequences of this interaction on the nonstationarity of the 
^ ' dephasing by considering a simple phenomenological model for a strongly coupled cluster of 
fluctuators. 

Let us thus consider a dissipative quantum two- level system (qubit) described by the 
Hamiltonian : 

n = ea'' + Aa^-^Xa\ (1) 
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where e and A arc control parameters which we assume to be time-independent and X de- 
scribes the fluctuating extra bias (noise) induced by the environment. In this letter, we focus 
on the case A = 0. However, extensions to the case A 7^ and a low frequency noise X may 
be treated along the lines of ref. [3]. Within a Dutta-Horn model [5,6], the noise X consists in 
the sum of contributions from independent fluctuators, modeled as random two-level systems 
coupled to an equilibrium bath {e.g. the phonons). This bosonic bath induces an elastic or 
electromagnetic coupling between the qubit and each fluctuator but also mutual couplings 0/ 
same origin between the fluctuators which are usually neglected. This interaction depends 
on the nature of the coupling to the bath : both elastic strain (phonons) and dipolar elec- 
tric coupling correspond to a long-range dipolar interaction [7, 8] The motivation to consider 
such couplings between the fluctuators comes in part from the analogous physics of TLS in 
amorphous solids: the relevance of these TLS couplings was stressed by Yu and Leggett [9], 
and found recent experimental evidence in the dielectric response of amorphous films [10]. 
These interactions are a possible source of broad distribution of relaxation times of the bath 
: in particular, using the language of pseudospins for the fluctuators, a strong ferromagnetic 
coupling between an ensemble of N spins forces these spins to tunnel simultaneously, there- 
fore reducing exponentially with the size N the tunneling amplitude of the cluster. Such a 
nonlinear relation between the size of a cluster and the corresponding relaxation time would 
naturally lead to a broad distribution for these times. Moreover, for randomly located fluc- 
tuators elastically coupled, we expect the interactions to be random and frustrating [7]. It 
is then natural to consider the noise produced by a mesoscopic ensemble of N interacting 
fluctuators [11]. We expect these interactions and the associated frustration to bring very 
slow relaxational dynamics for the cluster. However this problem is a formidable task, still 
unsolved in the case of classical spin glasses. We therefore adopt a phenomenological approach 
and focus on the expected general consequences of its physics. 

As in the case of independent fluctuators [5, 6] we will focus on a classical description of 
the noise X{t), considering an incoherent environment and neglecting back-action effects of 
the qubit. A common picture to describe the evolution of a glassy system is a random walk in 
phase space with broadly distributed time intervals r between the 'flips of the environment'. 
Hence, by analogy with the phase space trap models of classical glassy systems [13], we consider 
a general random intermittent noise, with randomly distributed heights Xi and durations 
of the plateaus, and waiting times between them (see Fig^l. For simplicity, we will focus 
on the case where the t] are much shorter than the r^, and given by their average rj. On time 
scales T ^ rj, the plateau is seen as a spike shifting the phase by XiT^. Moreover, to focus 
on dephasing and eliminate the drift of the accumulated phase of the qubit, the distribution 
V{x) is assumed to have zero mean and a finite width = (x^) (rj)^: It will be taken as 
Gaussian without loss of generality. Finally, we consider an algebraic distribution of waiting 
times Ti characterized by an exponent /i: 



where tq is a microscopic time. The slow algebraic decay of P{t) thereby implies divergences 
of all moments r" = dTP{T)T^'^ with n > ji. It turns out that three different classes of 
jjL have to be distinguished: (i) If both r and are finite [fx > 2), the dephasing induced 
by the noise in our model and within the usual model of Poissonian telegraph noise [5,6] are 
identical. Thus, we will call it the Poissonian class. However, in the case of slow dynamics, 
i.e. for /i < 2, the dephasing scenario differs considerably from the Poissonian class, (ii) For 
1 < /i < 2, when r is still finite, but already the second moment diverges, the waiting times 
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Ti start to fluctuate strongly around their average r and important corrections with respect 
to dcphasing due to Poissonian class arise, (iii) For < /i < 1, the divergence of the first 
moment t makes it even impossible to define a characteristic time between two consecutive 
spikes. In this case the two-point noise correlation function decays as 

x{t + t')x{t) - lot) X{t + t') ~ Y{if {t/t'f~^ , (3) 

where the overline denotes an average over realizations of X{t). Eq.Q implies that the power 
spectrum exhibits a divergence at low frequencies. Note that the explicit dependence 

of the correlator ((JJ on i implies a nonstationarity of the noise and consequently a time 
dependence of the amplitude of the noise [12]. In our model, a single source of noise with 
11 = 1 can lead to a l/f power spectrum whereas infinitely many independent fluctuators 
are necessary in the Dutta-Horn approach [5]. Within this latter approach, nonstationary of 
the noise comes from the choice of initial conditions for the slow fiuctuators. In our model, 
nonstationarity finds its origin in the very slow dynamics underlying the noise, not in a specific 
choice of an initial condition. 

X{t) -^0 

A I 

^ 

t' = 

Fig. 1 - Representation of a noise configuration X{t) 

The nonstationarity of the noise leads us to pay special attention to the dependence of 
dephasing on the age of the noise (denoted by tp below), i.e. a history dependence. To be 
definite we consider situations where the noise has been initialized at time t' = 0. The coupling 
to the qubit is turned on after a time tp, the preparation time of the qubit. Consequently, 
the two qubit states |±) accumulate a random relative phase between time tp and tp + t: 
$(ip,<) = J^''^* dt' X{t'). The corresponding dephasing factor is defined by the average over 
many realizations of noise 

i?(tp,t) = exp {i $(tp,t)) (4) 

and the dephasing time as its characteristic decay time: D{tp,T^) = e^^. In the following 
we will always consider preparation times tp that are large compared to the microscopic time 
tq. Under this assumption dephasing properties derived below do not depend on the imposed 
initial condition of the noise (at t' = 0). In particular, the noise is stationary for /x > 2, as in 
the case of usual Poissonian noise treated in [5,6]. 

Within the model J^J, the accumulated phase $(tp,t) performs a continuous time random 
walk (CTRW) [14]. The explicit dependence of D{tp,t) on tp can be accounted for using 
renewal theory [15]. The times Ti between two successive spikes being independent from each 
other (and thus of the history of the noise), the tp-dependence of ^{tp,t) and D{tp,t) are 
consequences of the tp-dependence of the distribution of ti , the time interval between tp and 
the first subsequent spike (see Fig. Indeed, the previous spike did not occur at tp but 
at some time tp — fi. Hence fi + ri is distributed according to P{t) in whereas ri 
is distributed according to a new distribution Ht^iri) which explicitly depends on tp. As 




4 



EUROPHYSICS LETTERS 



we show, this distribution Ht^{Ti) contains all information about the history of the noise. 
Following general ideas from renewal theory we separate noise configurations that have their 
first spike at tp + ri from the others to obtain an expression for Ht^{Ti) : 

HtM)^Pih + ^i)+ P(Ti+fi)5(ip-n)dfi, (5) 
Jo 

where S{t) denotes the renewal distribution, i.e. the density distribution of spike at time t. 
The behavior of S{t) follows from another renewal equation, 

S{t)^P{t)+ f P{T)S{t-T)dT, (6) 

Jo 

which states that a flip occurring at time t is either the first one, or that a previous flip 
occurred at time t — r, where t is distributed according to Denoting by L[f] the Laplace 
transform of the function /, Eq.lO can be rewritten as L[S] = L[P] / {1 — L[P]) . For jj, > 2, S 
is constant and Ht^{Ti) coincides with P{ti) at times tq. After a transcient regime, the 
noise has lost memory of its initial condition: the probability to flip to another state does not 
depend on the history of the noise. For fi < 2, the time dependence of S{t) implies an explicit 
history dependence of Ht^iji). For instance, for 1 < ^ < 2, S{t) ~ r~^(l + (ro/t)'*"^), and 
the average value of (ti) increases as t^~^. In the case /i, < 1 of diverging mean value, t — oo, 
the density of spikes S{t) decreases as and Ht^ differs therefore considerably from P. 

Now understanding the origin of aging, we decompose noise configurations into those 
without any event between tp and tp + i, and those with at least one event in the same interval, 
the first of which occurring at time t-p + t. In the latter case, the first event contributes a 

reinitializing simultaneously the noise at + t: 

D{tp,t) ^ lioih, t) + f{g) I Ht^{T)D{0,t-T)dT. (7) 

"'0 



In this equation, f{g) = ( e""^" ) = exp(— g /2) is the characteristic function of V{x) and 



Ho (ip , t) = Jf^°° dr Ht^ (t) denotes the probability that no spike occurs between tp and tp + t. 
More conveniently, denoting by L[Dt_p]{s) the Laplace transform of D{tp,t) considered as 
a function of t, and specializing ^ to tp = 0, we find L[Do]{s) = s~^(l — L[P]{s))/{l — 
f{g)L[P]{s)). Plugging it back into the Laplace transform of Q) we obtain 

This expression is one of the central results of this letter. Indeed, we can apply a numerical 
Laplace Transform inversion to obtain the complete behavior of the dephasing factor D{tp, t). 
Some of these results are shown in Fig. |21 The behavior of D{tp,t) is shown for /i = 0.8, 
H = 1.1 and both weak and strong coupling. 

We now turn to a discussion of the results we can derive exactly from Eq.®. As discussed 
above, the dephasing scenario differs qualitatively and quantitatively depending on the value 
of fj,. Generally, two regimes appear as a function of the coupling constant g, separated by a 
critical coupling constant gc'. While the dephasing time depends on g in the weak coupling 
regime, g < gc, it saturates for g > gc a.s a function of g (see Fig. EJ. As expected, it turns 
out that the critical coupling strength gc is of order 1 for the Poissonian class (/x > 2) as 
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Fig. 2 - Dephasing factor and obtained by numerical inversion of (HI for < 1 (left) and 1 < ^ < 2 
(right). For weak coupling, g < gdtp), the decay is exponential for t < r^. For times t > {/j. < 1) 
- t ^ {1 < fj, < 2) - and for strong coupling the decay is algebraic. For strong coupling and < 1 
the dephasing time exhibits a ip-dependence, which disappears as /i increases to higher values. Note 
the explicit tp-dependence of the critical coupling gc in the case ^ < 1. 



well as for 1 < /X < 2. However, in the case /j < 1, wo find that the range of the strong 
coupling regime increases with the age of the noise, i.e. gc decreases as a function of tp-. 
9c(tp) = '^(a')('''o/^p)''^^- This implies that any qubit surrounded by a noise with < /i < 1 
will eventually end up in the strong coupling regime. 

In the following we will discuss in detail the decay of D{tp,t) as a function of t and the 
scaling laws of as functions of g for the three different classes of /i and for both weak and 
strong coupling. In particular, we will compare our results with those for the Poissonian class. 
For jj, > 2, Ht (r) ~ P{t) and from © D{tp,t) ~ D{0,t). As expected, after a transcient 
regime {tp < tq), D{tp,t) becomes independant of tp. The results of [5,6] for a single fluctuator 
are then recovered for the Poissonian class. 

For 1 < fj, < 2, strong corrections with respect to the Poissonian class (/i > 2) occur in 
the strong coupling regime, g > 1. In that regime, the decay law of D{tp,t) is no longer an 
exponential but becomes algebraic D{tp,t) ~ [to/{to + t)]'^~^ ~ [to/(to + t + tp)]''"^. The 
dephasing time thereby exhibits an exphcit tp-dependence, ~ ro[l/e + (To/tp)^~^]~^/'^~^^ 
Note that this tp dependence is only important for values of [i close to one and disappears as 
[I approaches higher values (sec Fig.E)). For weak coupling, 5 < 1, the initial decay of D(tp, t) 
is very well described within a Gaussian approximation (second cumulant expansion). The 
decay is exponential, D{tp,t) = exp(— t/r^), and scales as ~ T^g^^ ■ Corrections with 
respect to Poissonian class are only subdominant for t < and disappear as fi increases to 
higher values. However, for times large compared to the dephasing time, t > log(r0/ro). 
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the above Gaussian approximation breaks down and the decay crosses over to a much slower 
power law (see Fig. O. 

In the case of < /.( < 1, the decay of coherence differs considerably from the Poissonian 
class in all regimes, due to the absence of a characteristic time scale in the waiting times 
distribution. For weak coupling g < gdtp) the decay for i < is accurately described by 
an exponential, D{tp,t) ~ exp[— (t/r^)^] with ~ roA(/i)g~^/^. Dephasing is therefore 
insensitive to the preparation time tp. This result can be recovered using the previously 
mcntionned Gaussian approximation. However, we emphasize the anomalous scaling of as 
a function of g (see Fig. But D{tp,t) is not fully described within this approximation. 
First of all, even at weak coupling, the Gaussian approximation breaks down for times t > t^. 
In this regime, the exponential decay is replaced by a much slower algebraic one: D(tp,t) ~ 
B{^) (t/r^)"^. This new behavior can be understood by noting that in this case and for 
t > D{tp,t) ~ D{tp,T^) no(r0,< — T^) ~ (t/r^— 1)"'^ . The second term of (j?)) corresponds 
to the anomalous random walk spreading of the phase it leads to the exponential decay (see 
above) which is subleading beyond r^. Hence, the leading term corresponds to the contribution 
of the noise configurations that did not change between and t. This situation of a main 
contribution induced by rare configurations is analogous to the physics of Griffiths singularities 
in disordered systems. In the strong coupling regime g > gdtp), D{tp,t) is dominated by the 
first term of Eq.jTI). In this limit, the distribution of the phase of the qubit starts spreading 
over [0, 27r]. Hence, most noise configurations produce a vanishing contribution to D{tp,t) 
and the whole average is dominated by those noise configurations that do not evolve during 
the experiment. The physics of this strong coupling regime is closely related to the physics 
of mean-field trap models of glassy materials [13], since the difference between quenched and 
annealed disorder is irrelevant in this case. The dephasing time then saturates as a function of 
g and becomes proportional to tp as seen in Fig. El Note that, at strong coupling, the Gaussian 
approximation breaks even down for short times and the decay of D(tp, t) is algebraic: it decays 
as D{tp,t) ~ 1 — sin(7r/x)/(7r^)(i/ip)^~^. For longer times, t > tp, it crosses over to a much 
slower decay, D{tp,t) ~ {tp/tY. 

Finally, in the marginal case /x = 1 of pure 1// noise the above analysis is confirmed quali- 
tatively. We mainly find logarithmic corrections to the above results, e.g. the critical coupling 

1/2 

scales as gc{tp) = [2(ro/ip) ln(ip/ro)] . For weak couphng, g < gdtp), the dephasing time 
is insensitive to tp and scales as ~ 2tq jln^j /g"^, whereas for strong coupling it depends 
algebraically on tp. 

In summary, we have studied decoherence of a qubit due to the noise generated by a 
slow collective environment. Within a simple phenomcnological model, we have derived exact 
expressions for the dephasing factor in various regimes. The crucial consequences of the non- 
stationarity of the noise are first the appearance of a history dependent coupling gdtp) which 
separate the weak and strong coupling regimes and, second, a non exponential decrease of the 
dephasing factor D{tp,t). Depending on the broadness (fi) of the distribution of relaxation 
times of the environment, this nonexponential decay appears either before {fi < 1) or after r^, 
but is a clear signature of the nonstationarity of the noise {i.e. /x < 2 in our model). In a more 
refined description of a collective noise source consisting of a collection of coupled clusters, 
we expect the nonstationarity of the dephasing to be related to the number of contributing 
clusters. Similarly to the studies of Paladino et al. [5] for usual telegraphic fluctuators, the 
nonstationary dephasing in our model will survive provided the phase $(ip, t) is dominated by 
a few strongly coupled slow clusters. In the opposite limit of many contributing clusters, the 
usual Gaussian stationary dephasing should be recovered. Thus the search for such nonsta- 
tionarity in the dephasing of simple solid state qubits would be of main interest for the correct 
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characterization of the source of 1// noise in these samples. Note however that dephasing is 
experimentaUy studied by repeated interference experiments without noise reinitiahzation in 
between. The dephasing is then characterized by a time-average instead of the average over 
configurations D(tp,t). In the case of a nonstationary noise these two expressions wiU differ 
(non ergodicity for < /i < 1) and special care should be given in analyzing the experimental 
results. 

J. Schriefl thanks Yu. Makhlin for very useful discussions. P. Dcgiovanni thanks the Insti- 
tute for Quantum Computing (Waterloo) and Boston University for support and hospitality 
during completion of this work. 
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